mbSet<-Init.mbSetObj()
mbSet<-SetModuleType(mbSet, "mdp")
mbSet<-ReadSampleTable(mbSet, "metadata 6-Yes and No.txt");
mbSet<-Read16STaxaTable(mbSet, "taxonomy_19.01.txt");
mbSet<-Read16SAbundData(mbSet, "ASV-feature-table_19.01.txt","text","Others/Not_specific","T","false");
mbSet<-SanityCheckData(mbSet, "text","sample","true");
mbSet<-SanityCheckSampleData(mbSet);
mbSet<-SetMetaAttributes(mbSet, "1")
mbSet<-PlotLibSizeView(mbSet, "norm_libsizes_0","png");
mbSet<-SanityCheckData(mbSet, "text","none","false");
mbSet<-SanityCheckSampleData(mbSet);
mbSet<-SetMetaAttributes(mbSet, "1")
mbSet<-PlotLibSizeView(mbSet, "norm_libsizes_1","png");
mbSet<-CreatePhyloseqObj(mbSet, "text","Others/Not_specific","F" , "false")
smpl.nm.vec <- c("RS001.AT","RS002.AT","RS003.AT","RS005.AT","RS007.AT","RS008.AT","RS009.AT","RS010.AT","RS011.AT","RS012.AT","RS013.AT","RS014.AT","RS015.AT","RS016.AT","RS017.AT","RS019.AT","RS020.AT","RS021.AT","RS022.AT","RS023.AT","RS024.AT","RS025.AT","RS026.AT","RS027.AT","RS028.AT","RS029.AT","RS030.AT","RS031.AT","RS032.AT","RS033.AT","RS035.AT","RS036.AT","RS037.AT","RS038.AT","RS039.AT","RS040.AT","RS042.AT","RS043.AT","RS044.AT","RS045.AT","RS046.AT","RS048.AT","RS049.AT","RS050.AT","RS051.AT","RS052.AT","RS053.AT")
mbSet<-UpdateSampleItems(mbSet);
mbSet<-ApplyAbundanceFilter(mbSet, "prevalence", 0, 0.1);
mbSet<-ApplyVarianceFilter(mbSet, "iqr", 0.0);
mbSet<-GetLibscale(mbSet);
mbSet<-PerformNormalization(mbSet, "none", "none", "none", "true",1427);
mbSet<-PerformBetaDiversity(mbSet, "beta_diver_0","PCoA","bray","expfac","type","none","OTU","","Chao1", "yes", "adonis", "png", 72, "default", "true");
mbSet<-PCoA3D.Anal(mbSet, "PCoA","bray","OTU","expfac","type","","Chao1","beta_diver3d_0.json")
mbSet<-PerformBetaDiversity(mbSet, "beta_diver_1","PCoA","bray","expfac","type","none","OTU","","Chao1", "yes", "adonis", "png", 72, "default", "true");
mbSet<-PerformCategoryComp(mbSet, "OTU", "adonis","bray","type","true");
mbSet<-PCoA3D.Anal(mbSet, "PCoA","bray","OTU","expfac","type","","Chao1","beta_diver3d_1.json")
mbSet<-PerformCategoryComp(mbSet, "OTU", "adonis","bray","type","true");
mbSet<-PlotAlphaData(mbSet, "filt","alpha_diver_0","Chao1","type","OTU", "default", "png");
mbSet<-PlotAlphaBoxData(mbSet, "alpha_diverbox_4","Shannon","type","default", "png");
mbSet<-PerformAlphaDiversityComp(mbSet, "nonpar","type","true");
mbSet<-PlotAlphaData(mbSet, "orig","alpha_diver_5","Observed","type","OTU", "default", "png");
mbSet<-PlotAlphaBoxData(mbSet, "alpha_diverbox_5","Observed","type","default", "png");
mbSet<-PerformAlphaDiversityComp(mbSet, "nonpar","type","true");
mbSet<-PlotAlphaData(mbSet, "orig","alpha_diver_6","Chao1","type","OTU", "default", "png");
mbSet<-PlotAlphaBoxData(mbSet, "alpha_diverbox_6","Chao1","type","default", "png");
mbSet<-PerformAlphaDiversityComp(mbSet, "nonpar","type","true");
mbSet<-PlotTaxaAundanceBar(mbSet, "taxa_alpha_0","Phylum","type", "null", "barraw",10, "set3","sum",10, "bottom", "F", "png");
mbSet<-PlotTaxaAbundanceBarSamGrp(mbSet, "taxa_alpha_1","Phylum","type", "none", "barraw",10,"set3","sum", 10, "top", "F", "png");
mbSet<-PlotTaxaAbundanceBarSamGrp(mbSet, "taxa_alpha_2","Phylum","type", "none", "barnorm",10,"set3","sum", 10, "top", "F", "png");
mbSet<-PrepareHeatTreePlot(mbSet, "type", "Genus", "dbgr", "dft", "No_vs_Yes", 0.05, "false", "all", "true", "heat_tree_0","png");
mbSet<-PerformRNAseqDE(mbSet, "EdgeR",0.05,"type","NA","OTU",0.0, "null1", "null2");
mbSet<-PerformRNAseqDE(mbSet, "EdgeR",0.05,"type","NA","OTU",0.0, "null1", "null2");
mbSet<-GenerateCompJson(mbSet, "comp_graphics_0.json","png" ,"EdgeR",1,"OTU", "Phylum",0.05);
mbSet<-GenerateCompJson(mbSet, "comp_graphics_1.json","png" ,"EdgeR",1,"OTU", "Phylum",0.05);
mbSet<-PerformRNAseqDE(mbSet, "EdgeR",0.05,"type","NA","Genus",0.0, "Yes", "No");
mbSet<-GenerateCompJson(mbSet, "comp_graphics_2.json","png" ,"EdgeR",1,"Genus", "Phylum",0.05);
mbSet<-PlotBoxData(mbSet, "box_plot_0","g__Bacillus","boxplot", "png");
mbSet<-PerformRNAseqDE(mbSet, "EdgeR",0.05,"type","NA","Species",0.0, "Yes", "No");
mbSet<-GenerateCompJson(mbSet, "comp_graphics_3.json","png" ,"EdgeR",1,"Species", "Phylum",0.05);
mbSet<-PlotBoxData(mbSet, "box_plot_1","Lactobacillus_johnsonii","boxplot", "png");
